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Abstract 

Residual late-time dark matter particle annihilations during and after Big-Bang Nucleosyn- 
thesis (BBN) may alter the predicted cosmological abundances of the light elements. Within the 
constrained minimal supersymmetric extension of the Standard Model (the CMSSM) with a neu- 
tralino LSP, we find negligible effects on the abundances of Deuterium, 3 He, 4 He and Li predicted 
by homogeneous BBN, but potentially a large enhancement in the predicted abundance of 6 Li. 
This enhancement may be as much as two orders of magnitude in the focus-point WMAP strip 
and in the coannihilation and funnel regions for large tan/5 for small Tn\/2, and the effect is still 
significant at large "2-1/2- However, the potential 6 Li enhancement is negligible in the part of the 
coannihilation strip for tan /3 = 10 that survives the latest LHC constraints. A similar enhancement 
of the 6 Li abundance may also be found in a model with common, non-universal Higgs masses (the 
NUHM1). 
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1 Introduction 



The success of homogeneous Big-Bang Nucleosynthesis (BBN) is one of the lynchpins of 
cosmology. Using the baryon-to-photon ratio, 77, inferred from measurements of the cosmic 
microwave background (CMB) radiation [T], homogeneous BBN predicts successfully the 
astrophysical abundances of Deuterium, 3 He and 4 He [2H7]. On the other hand, there are 
issues with the abundances of 7 Li [7J and potentially with 6 Li [8]. In particular, the predicted 
abundance of 7 Li is considerably larger than the range suggested by observations [9HT8] 0, 
and there are suggestions that the astrophysical value of the 6 Li abundance may be much 
higher than predicted by homogeneous BBN [T3]. However, one should note that the line 
asymmetries which have been interpreted as 6 Li could be the result of convective processes 
affecting 7 Li [20]. 

We and others have investigated previously whether the late decays of massive particles, 
such as the gravitino in the constrained minimal supersymmetric extension of the Standard 
Model (the CMSSM) with a neutralino as the lightest supersymmetric particle (LSP), could 
improve significantly the 6 Li and 7 Li abundances predicted by homogeneous BBN [2TH45]. 
We did not find a solution to the 6 Li problem, but we did find a region of supersymmetric 
parameter space where gravitino decays might alleviate or even solve the 7 Li problem [I5H6] . 
On the other hand the 7 Li problem might have a more banal solution, such as the existence 
of a suitable carbon, boron or beryllium resonance [17] . 

In addition to decays, the late-time annihilations of cold dark matter may also affect 
the abundances of the light elements [25j [35l EH1 SHI SS] - I n particular, these annihilations 
may have a significant effect on the abundance of 6 Li [36]. There it was argued that 6 Li 
production may occur if the s-wave annihilation cross-section is sufficiently large, and it 
was assumed that the relic density of the annihilating dark matter particles is controlled 
largely by the s-wave part of the cross-section. However, in supersymmetric models where 
the LSP is a neutralino, such as the CMSSM, the relic density is in fact largely determined 
by the p-wave part of the cross-section, which by the time of BBN is essentially ineffective. 
Therefore, a re-analysis of the suggestion of [36] in the context of the CMSSM and related 
models is timely, and is the subject of this paper. 

In this paper we study the possible effects on the cosmological light-element abundances 
of residual late-time annihilations of neutralino LSPs during or after BBN [35 | l36 | H8] Fl We 

1 A globular cluster star with a 7 Li abundance comparable to the BBN prediction has recently been 
observed [19] : this value may be due to production by a previous generation of stars. 

2 Recent papers have also considered the BBN consequences of WIMP models having residual annihilations 
increased by Sommerfcld or Briet-Wigncr enhancements [50]. In the CMSSM model, neither effect occurs, 
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find negligible effects on the abundances of Deuterium, 3 He, 4 He and 7 Li predicted by ho- 
mogeneous BBN, but potentially a large enhancement in the predicted abundance of 6 Li, as 
suggested in [36]. The physics of this effect is the following. It is well understood that the fa- 
mous A = 5 gap in the spectrum of stable nuclei impedes the production of heavier nuclei in 
BBN. The dominant mechanism for making 6 Li in annihilating-particle scenarios is initiated 
by p and n spallation of 4 He. This yields many A = 3 nuclei with only a tiny reduction in 
4 He abundance. The tritium and 3 He nuclei are produced with large, nonthermal energies, 
and subsequently slow down due to ionization losses, but have some probability of inducing 
t(a,n) 6 Li or 3 He(a,p) 6 Li reactions first. In this way, an amount of 6 Li may be produced 
that is large relative to the standard homogeneous BBN abundance, without making large 
amounts of extra Deuterium and A = 3 or reducing the 4 He abundance, and leaving the 7 Li 
abundance unaffected. 

In the CMSSM [5TJ, it is assumed that all supersymmetry-breaking gaugino masses have 
a common value mu2 at some grand unification scale before renormalization, and likewise 
all the soft supersymmetry-breaking scalar masses are assumed to have a common m . The 
other parameters of this model are the (supposedly) universal trilinear parameter Aq (taken 
here to be Aq = 0), and the ratio of Higgs vacuum expectation values, tan/3. In addition, one 
must specify the sign of the Higgs mixing mass, /i, which is generally taken to be positive 
in the CMSSM so as to improve compatibility with measurements of — 2 and 6—7-37 
decays. As is well known, in the CMSSM there are strips in the (m 1( / 2 ,mo) planes for fixed 
Aq and tan j3 along which the relic density of the neutralino LSP, \i nes within the range 
favoured by WMAP and other astrophysical observations [52] . At relatively low values of mo 
there is generically a coannihilation strip close the boundary where the LSP would become 
charged that extends, for large tan/3, into a funnel at large mi/2 where the LSPs annihilate 
rapidly via direct-channel heavy Higgs resonances. At relatively high values of m , there 
is a focus-point strip close to the boundary of consistent electroweak symmetry breaking. 
In light of present experimental constraints from the LHC and elsewhere, plausible values 
of tan/3 range between ~ 10 and ~ 55 [53]. Representative (mi/2,^0) planes for these 
values of tan /3 are shown in Fig. [TJ and are discussed below in more detail. In this paper, 
we explore the effects on cosmological light-element abundances of residual late-time XX 
annihilations for CMSSM parameters along the coannihilation/funnel and focus-point strips 
for these reference values of tan /3, and also remark on additional possibilities in one- and two- 
parameter generalizations of the CMSSM with non-universal soft supersymmetry-breaking 
contributions to Higgs masses (the NUHM1 and NUHM2) [55ti57]. 

due to the lack of a light boson and of extreme degeneracy in the funnel, respectively. 
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Figure 1: Left: The CMSSM (m^jmo) plane for A = and tan/3 = 10, and Right: 
the corresponding plane for tan /3 = 55, both with \x > 0, displaying contours of the 6 Li 
abundance including the effects of late-time XX annihilations. Contours of the 6 Li abundance 
are coloured light blue, and the WMAP- compatible strips of parameter space are shaded 
dark blue. The brown shaded region at large and small mo is excluded because the LSP 
would be charged, and in the pink shaded region at small m\/2 and large m there would be 
no consistent electroweak vacuum. Also shown are the exclusion by LEP searches for the 
Higgs boson (red dash-dotted line) and charginos (black dashed line), and by LHC searches 
for sparticles (purple solid line) /53[[5ffi. The green shaded region is excluded byb^- sj, and 
the paler pink region is favoured by g^ — 1 at the 1- (2-)o~ level, as indicated by the dashed 
(solid) black lines. 

We find that such late annihilations have no significant effect on the cosmological abun- 
dances of Deuterium, 3 He, 4 He and 7 Li in any of the CMSSM scenarios studied. However, 
they may enhance the 6 Li abundance by up to two orders of magnitude: in some instances, 
we find 6 Li/H ~ 10~ 12 , compared to the value ~ 10~ 14 found in standard homogeneous 
BBN [58,59]. This possibility arises at low m\ji along the focus-point strips for tan/3 = 10 
and 55, and also along the coannihilation strip for tan/3 = 55. The values of 6 Li/H along 
these strips decrease to ~ 10~ 13 at large min, which is typical also of values in the funnel 
region for tan/3 = 55. On the other hand, we find no substantial enhancement of 6 Li/H 
along the coannihilation strip for tan/3 = 10 for values of m 1( / 2 consistent with the LHC 
constraints [531 153], though values as large as ~ 10~ 13 might have been reached at lower 
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mi/2. The possibility of a large 6 Li enhancement is extended in the NUHM1 to a large range 
of mo with tan/3 = 10 and low mi/2- 

2 The Cosmological Lithium Problems 

Cosmologically, the dominant Lithium isotope is 7 Li, whose abundance is commonly inferred 
from observations of low-metallicity halo dwarf stars. These indicate a plateau of Lithium 
versus metallicity [9], with 

= (1-23123) x ID" 10 , (1) 

n / halo* 

whereas observations of globular clusters [12] yield somewhat higher values: 



7 Lf 

— ) = (2.35 ±0.05) x 10" 10 . (2) 
H 'GC 



For comparison, the standard BBN result for 7 Li/H is (5.12l ; 62 ) X 10~ 10 [TJ, and the difference 
between this and (JU [2]) constitutes the cosmological 7 Li problem. As already mentioned, this 
might be resolved by new physics beyond the Standard Model such as late decays of massive 
gravitinos jl6], or by some undocumented Standard Model effect such as a suitable carbon, 
boron or beryllium resonance [47J. The cosmological 7 Li problem is not our focus in this 
paper. 

6 Li has been observed in some halo stars [60] with [Fe/H] ~ -2, and with an isotopic ratio 
that is 

(^) ~ 0.05. (3) 

\ / halo* 

These observations are consistent with the results of Galactic cosmic-ray (GCR) nucleosyn- 
thesis [nnKnUEZj , though see below for results at lower metallicity. This confirms that most 
of the Lithium is in the form of 7 Li, leaving unscathed the cosmological 7 Li problem. 

However, a recent paper has reported the presence of a similar isotopic abundance in halo 
stars over a broad range of metallicities that extends to significantly lower values ([Fe/H] ~ -1 
to -3) [TI1I60] . The inferred plateau 6 Li/H ratio ~ (6 to 25) x 10~ 12 is about 1000 times higher 
than the 6 Li/H ratio predicted by standard homogeneous BBN [58j[59], namely 6 Li/H ~ 
10~ 14 . The isotopic ratio (131) cannot be explained by conventional GCR nucleosynthesis, 
at the lowest metallicities: this is the cosmological 6 Li problem. The reliability of the 6 Li 
plateau at very low metallicity has been questioned [20], so the 6 Li problem should be taken 
with a grain of salt. But in any case, these exiting if controversial results demonstrate that 



4 



6 Li abundances at levels 6 Li/H < few x 10~ 12 are at or near the reach of present observational 
techniques. 

Thus the current observational situation is evolving, but without question is interesting: 
at the very least, the present results serve as upper limits to primordial 6 Li, and impose 
bounds on nonstandard BBN. At most, current data may already point to a primordial 6 Li 
problem which would demand new BBN physics, and probe its details. Our focus in this 
paper is to determine the 6 Li production and its observational implications in the context of 
some of the most popular supersymmetric dark matter scenarios. 

It has been proposed that some decaying-particle scenario might produce 6 Li at the 
plateau level with some destruction of 7 Li [H1I63H6T] . offering the possibility of solving both 
Lithium problems simultaneously. However, we note that solving the 6 Li problem would 
use up only a small fraction of the 7 Li whose destruction would be needed to solve the 7 Li 
problem, leading one to consider separate solutions for the two Lithium problems. It is 
also possible that the 6 Li problem might be explained by nucleosynthesis due to cosmolog- 
ical cosmic rays produced at the epoch of structure formation [8|IS8]. We have previously 
demonstrated that late-decaying massive gravitinos might resolve the 7 Li problem within 
the CMSSM framework [16]. Here we show that the 6 Li problem might, independently and 
in parallel, have at least a partial supersymmetric solution, via the late annihilations of 
neutralino LSPs. 



3 Residual Late-Time Neutralino Annihilations 

Assuming that the lightest neutralino x is the LSP, and that R-parity is conserved, the 
relic neutralino density is essentially fixed at a freeze-out temperature Tf ~ m x /20. At 
lower temperatures, the local density of neutralinos, n x , decreases as the universe expands 
(presumably) adiabatically, and subsequent annihilations have very little effect on the dark 
matter density, but may have important effects on the light-element abundances 



The rate per volume of annihilation events is 

<?ann = Mann, ( 4 ) 

and so the annihilation event rate per x is 

Tann = — = ^(<™)ann% = \ (<™)ann^ X n b, (5) 
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and thus the annihilation event rate per baryon is 

= r ann F x = ^(o-w) ann F x 2 n b , (6) 

Tib ^ 

where the x abundance is 



rib 



At the temperatures of interest here, Tbbn ^ 1 MeV <C m x , the annihilation rate coefficient 
(c rt, )ann is very well approximated as a constant, the value of which depends on the specific 
underlying supersymmetry model. In (J3J), we are interested in only the s-wave part of the 
cross-section whereas a combination of s- and (mainly) p-wave cross-sections is constrained 
by the requirement of reproducing the present dark matter density within errors. 

The annihilations inject nonthermal Standard Model particles, including both electro- 
magnetic as well as hadronic species. For electromagnetic products we need only track the 
total energy injected per annihilation. For nonthermal hadrons (nucleons) h = n,p, we calcu- 
late the spectrum Q^ nn (e) of annihilation products, normalized such that J Qf^ n (e) de = B^, 
the expected number of h created per annihilation. Then the injection/source rate of h due 
to annihilations, per unit volume, per unit time, and per unit kinetic energy e, is 

j iw-ann 

hM ~ q,nnQr(e) . (8) 



dV dt de 



These particles then lose energy as they propagate in the cosmic plasma. The propagated 
spectrum of nonthermal particles must be calculated, and this produces the nonthermal 
reactions on ambient thermal light nuclides that perturb BBN. 

The effect of nonthermal particle injection in BBN has been well-studied in the case of 
decays of some unstable particle X. Much of the physics carries over here, once one makes 
the appropriate substitution of abundances nx — > n x and of annihilation rate for decay rate: 
Tx = t^ 1 — > r ann . After injection, the nonthermal particle propagation remains the same, 
and we treat this as in |45j . We also adopt the same set of nonthermal BBN reactions which 
we include in the same manner, making the appropriate substitution of annihilations for 
decays. 



3.1 Order-of-Magnitude Calculation 

Before turning to our numerical results, we first present an order-of-magnitude calculation 
that illustrates the basic physics in play, and also serves as a check on our numerical results. 
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The total number of annihilation events per baryon occurring after a given time U is the 
time integral of eq. (jSJ) 

A4nn = / r ann Y x dt ~ Y^(av) ann n h (ti) U (9) 



( Mann \ /300GeVY 
V10" 26 cm 3 s~V V m x ) 



5 x 10 events/baryon . (10) 



Our fiducial values correspond to tj ~ 100 sec and T, ~ 100 keV, since this marks the epoch 
when the 4 He abundance becomes large. 

Given this number of annihilations per baryon, we now need the branching for 6 Li pro- 
duction per annihilation. As discussed earlier, nonthermal particles from annihilations or 
decays produce 6 Li as a secondary by-product of 4 He spallation: 



^nonthermal He Y H non ^ nerma ] 

1 nonthermal 



3 H non thermal + 4 He 6 Li + U , (11) 



and similarly with nonthermal 3 He; nonthermal D also contributes but is subdominant. As 
discussed for the late-decay case in [15], each late annihilation produces a mass-3 abundance 
increment AY( 3 A) which is given in the thin target limit by 

a(i\a — > inelastic) 

where Bn ~ 0.4 is the number of nucleons per annihilation. Typically this increases the mass- 
3 abundance by an amount AY( 3 A) ~ 10 -9 1bbn( 3 ^4) ~ 10~ 5 , i.e., much smaller than 
the standard primordial abundance, and thus we do not expect substantial perturbations to 
mass-3 nuclides, or to D, which has similar cross sections, or to 4 He. 

The energetic A = 3 particles are slowed in the cosmic plasma by ionization and related 
losses, with a range i?3 = J(dE/dX)~ 1 dE, where dE/dX is the loss rate per thickness 
dX = p h dx in [g/cm 2 ]. Hence the stopping length is Rz/ Ph- The fraction of mass-3 nuclides 
which produce 6 Li before stopping is this stopping length divided by the mean free path for 
6 Li production, namely: 

f{ 3 A -> 6 Li) ~ n a a( 3 Aa 6 Li)^ ~ y ^Aa^^R, ^ ? ^ 1Q _ 4 

Ph m h 

Collecting these results, the residual late-time annihilation contribution to the 6 Li abundance 
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per baryon is 



Ay( 6 Li) = AY( 3 A) f( z A -> 6 Li) (14) 
„ ,,2/ v °( N ® ->• 3A + ■ • • ) cr^Aa -»■ 6 Li)i? 3 

~ 5 Jv r Y 2 ( C T W ) ann F a 4- — r-rr — — «b(*i) *i (15) 



The numerical results given above are evaluated for ^ = 100 sec, and we also take 
i?3 = lg/cm 2 and a( 3 Aa — » 6 Li) = 30 mb. This formula gives the scaling Ay( 6 Li) oc 
B^(av) fmri {VL x /m x ) 2 which we verify with our full numerical results, and the normaliza- 
tion agrees to within a factor ~ 2. This agreement lends confidence in our code and our 
understanding of the physics. 

3.2 Numerical Results 

We turn now to our full numerical results. In order to establish the context for our subsequent 
analysis of the possible annihilation effects along the strips in CMSSM parameter space that 
are compatible with WMAP and other constraints on the present-day dark matter density, 
we first discuss the full CMSSM (mi/2, mo) planes shown in Fig. [TJ The light blue lines 
are contours of the 6 Li abundance, and the relic density is WMAP-compatible pQ along the 
dark blue strips, assuming that the lightest neutralino \ is t ne LSP and is stable, as in 
.R-conserving models. There would be no consistent electroweak vacuum in the pink shaded 



SP in the brown shaded region, 
Regions to the left of the red 



region at small mi/2 and large mo, the lighter f would be the 
and the green shaded region is excluded by b — > s-f decay 
dash-dotted (black dashed) (purple) line are excluded by searches for the Higgs boson at LEP 
(charginos) (LHC searches for sparticles) [531154] . In the paler pink region the supersymmetric 
contribution remedies the discrepancy between the experimental measurement of — 2 and 
theoretical calculation within the Standard Model using low-energy e + e~ data, with 1- (2-)cr 
consistency being indicated by the dashed (solid) black lines. 

We see in the left panel of Fig. [T] showing the (m.1/2, mo) plane for tan (3 = 10 that most 
of the lower (coannihilation) WMAP strip has 6 Li/H < 10 -13 , whereas the upper (focus- 
point) strip may have 6 Li/H as large as 10 -12 . There is a region where 6 Li/H seems able to 
exceed 10" 11 , but this is well inside the region between the WMAP strips, where the relic 



3 According to conventional Big-Bang cosmology and in the absence of R violation, the LSP x would be 
ovcrdcnsc in the regions between the WMAP strips. It would be undcrdense in the regions between these 
strips and the pink and brown shaded regions. 



s 



X is overdense according to conventional Big-Bang cosmology. In the right panel of Fig. [T] 
for tan/3 = 55, we see that 6 Li/H ~ 10~ 13 along the coannihilation strip and in the funnel 
region at large mm and mo where the relic density is brought into the WMAP range by 
rapid annihilations through direct-channel H/A resonances, though somewhat larger values 
of 6 Li/H are possible at small mm. Along the focus-point strip, we see that values of 6 Li/H 
~ 10~ 12 are also possible at small mm, falling to ~ 10~ 13 at large mm. The range 6 Li/H 
~ 10~ n is never attained for tan/3 = 55, even in the overdense region of the (mi/2, wio) 
plane. 

We now focus on the WMAP strips in Fig. [U The left panel of Fig. [2] displays the figure- 
of-merit combination (av) ann (Q x h 2 / m x ) 2 as a function of mi/2 along the WMAP strips in 
the CMSSM for tan /3 — 10 and 55. We see that (av) ann (Q x h 2 /m x ) 2 along the coannihilation 
strip for tan /3 = 10 is much smaller than along the other strips. This can be understood from 
the fact that along this strip several coannihilation processes involving sleptons contribute to 
reducing the relic \ density into the WMAP range, and that their relative contributions be- 
come more important as mm increases. In addition, along this strip the s-wave cross-section 
relevant during BBN is significantly smaller than the p-wave cross-section that dominates 
during freeze-out. These coannihilation processes are less important along the corresponding 
strip for tan/3 = 55, and unimportant in the funnel region and along the focus-point strips, 
where the s-wave cross-section becomes comparable to the total cross section. Hence, along 
these strips (crv) aan must be larger, in order to bring the relic density down into the WMAP 
range unaided. 

The XX annihilations feed many different particle species into the cosmological back- 
ground, initially with nonthermal spectra that we model using PYTHIA [69]. The only species 
that survive long enough to interact significantly with background nuclei are protons and 
neutrons (and their antiparticles) and photons. The former are far more important for the 
nuclear reactions of interest here, so we focus on their numbers and spectra. The right 
panel of Fig. |2] displays the numbers of protons (solid or dotted lines) and neutrons (dashed 
or dash-dotted lines) produced per annihilation event, again along the WMAP strips for 
tan /3 = 10, 55 discussed previously. We see that in general the numbers of protons and neu- 
trons increase significantly as mm increases, with some bumps as new annihilation thresholds 
are crossed. 

Fig. 3 of [70] displays the most important branching fractions for final states in XX 
annihilations as functions of mm along the WMAP strips for tan/3 = 10 and 55, which 
include the final states r + r~, bb, W + W~, ti, hZ and ZZ. Of these, the r + r~ final state 
clearly yields no baryons, while the numbers of baryons yielded by the final states W + W~, 
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Figure 2: Left: The figure of merit for the late-time annihilation rate, namely the 
velocity-averaged XX annihilation cross section divided by the square of the neutralino mass, 
(cn;) ann (f2 x /i 2 /m x ) 2 , along the WMAP strips in the coannihilation, focus-point and funnel 
regions for tan f3 = 10, 55, Aq = and \x > 0, as functions of m\/2- We see that the figure of 
merit along the tan (3 = 10 coannihilation strip is much smaller than along the other strips, 
and that all decrease rapidly as m^j increases. Right: The numbers B PjU of protons (solid or 
dotted lines) and neutrons (dashed or dash-dotted lines) produced per XX annihilation event, 
as calculated using PYTHIA, along the WMAP strips in the coannihilation, focus-point and 
funnel regions for tan/3 = 10,55, Aq = and \i > 0, as functions o/m^. We see that in 
general the numbers increase significantly as mi/2 increases. 



hZ and ZZ are all independent of the annihilation centre-of-mass energy 2m x . Only the bb 
and ti final states yield numbers of baryons that increase with the annihilation centre-of-mass 
energy 

Fig. [3] displays the spectra of protons (upper panel) and neutrons (lower panel) for the 
W + W~, hZ and ZZ final states for m x = 250 GeV, and for the bb final state for m x = 
100, 250 GeV, all calculated using PYTHIA. We display the number of protons or neutrons 
per unit of the parameter x = \/ Ef — mf/m XJ where % = p,n. The proton and neutron 
spectra are almost identical. They differ in the small x region primarily because of the 
difference of m p and m n . We also see that the W + W~, hZ and ZZ final states yield rather 
similar spectra, with the spectrum from the hZ final state rising slightly higher. The spectra 
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of baryons from bb final states rise from being lower at m x = 100 GeV to being higher at 
m x = 250 GeV. 

Together with Fig. 3 of [70], Fig. |3] enables us to understand the salient features of the 
baryon production rates shown in the right panel of Fig. |2j The large branching fraction for 
t + t~ suppresses baryon injection along the coannihilation strip for tan j3 = 10, particularly 
for small m^j but less so for large mx/2 where the W + W~ branching fraction grows. 

Following their injection into the primordial plasma, some of the nucleons cause spallation 
of 4 He, yielding A = 3 nuclei as discussed above. These are produced with large, nonthermal 
energies and subsequently thermalize, but may previously induce t(a,n) 6 Li or 3 He(a,p) 6 Li 
reactions. Fig. H] displays the enhancement of the cosmological 6 Li abundance that we find 
along the WMAP strips discussed above. The homogeneous BBN value ~ 10~ 14 is attained 
at large m 1( / 2 along the WMAP coannihilation strip for tan/3 = 10, but much larger values 
are possible along the other WMAP strips ^j, where we find 6 Li/H ~ 10~ 13 at large mx/2 to 
10~ 12 at small mx/2. We have found that the enhancement of 6 Li scales very closely with the 
combination B N {crv) a _ nn /m x , as was to be expected. 

We recall that the enhancement of 6 Li/H that would be required for consistency with (j3J) 
is by a factor ~ 1000, rather than the factor of up to ~ 100 that we find here. However, as 
we have noted there remains a question as to whether or not the plateau ratio of 0.05 should 
be attributed to 6 Li. The abundance of 6 Li we find here is potentially observable and would 
in fact be seen as a plateau extending to low metallicities. Optimistically, we could envision 
6 Li observations playing a role in discerning between supersymmetric models. In any case, 
we regard the enhancement we find as already an interesting contribution to the analysis of 
the 6 Li problem. 

3.3 Exploration of Non-Universal Higgs Models 

It is quite possible that some modifications of the CMSSM might yield even greater enhance- 
ments of the 6 Li abundance. To be successful in this respect, it is apparent from (fl6l) that 
such a model would require a relatively large annihilation cross section (av) ann combined 
with a small value of m x , as in the focus-point region of the CMSSM. There, the relatively 
large value of (crt>) ann is made possible by the admixture of a Higgsino component in the x, 
and along this strip the low value of m x is consistent with the LHC and other constraints [53J. 

In an initial probe of other possibilities for a large enhancement of the 6 Li abundance, we 
have explored the NUHM1 model, in which the soft supersymmetry-breaking contributions 



4 The region of enhanced 6 Li along the WMAP coannihilation strip for tan j3 = 10 with m-i/2 < 400 GeV 
is now excluded by the unsuccessful LHC searches for supersymmetry [53 1 154], 
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Figure 3: The spectra of protons (upper panel) and neutrons (lower panel) injected by XX 
annihilations into the Zh, W + W~ , ZZ and bb (for m x =100 and 250 GeV) final states, as 
calculated using PYTHIA. 
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Figure 4: The enhancement of the cosmological 6 Li abundance as a function of mi/ 2 along 
the WMAP strips discussed in the text. The standard homogeneous BBN value ~ 10 -14 is 
attained at large mi/2 along the WMAP coannihilation strip for tan/3 = 10. 

to the Higgs masses have a common value that differs from mo- It is known that in this 
model the Higgsino component in the LSP \ may be enhanced at values of and m 
away from the focus-point region, thanks to a level-crossing transition at particular values 
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of /i/mi/2 [56]. In the CMSSM, the value of /i is generally fixed by applying the conditions 
for a consistent electroweak vacuum. However, in the NUHM1 the value of /V m i/2 can 
be adjusted by varying the degree of non-universality in the soft supersymmetry-breaking 
Higgs masses, enabling a WMAP-compatible relic density to be found in models with values 
of (mi/2, 7n ) different from those allowed in the CMSSM. 

We have explored the conditions under which such transition regions in the NUHM1 
may yield an enhancement of 6 Li/H similar to, or (possibly) greater than the value ~ 10~ 12 
attainable in the CMSSM in the focus-point region. To this end, we have studied over a dozen 
NUHM1 parameter planes. In no case did we find enhancements of 6 Li/H significantly larger 
than in the CMSSM (and this is also the case in some planes we explored in the NUHM2, 
in which both Higgs soft supersymmetry-breaking masses are treated as free, non-universal 
parameters) . 

Fig. [5] shows results in a couple of selected NUHM1 parameter planes. The left panel 
shows an (mi/2,777.0) plane for tan/3 = 10, Aq = and fixed \x = 250 GeV. In this case, 
there is a near- vertical WMAP-compatible strip in a transition region at ~ 400 GeV. 
This transition strip is compatible with the LEP Higgs constraint, and the upper part of the 
strip above mo ~ 700 GeV is compatible with the constraints imposed by LHC searches for 
sparticles. We see that the 6 Li/H ratio is remarkably constant at ~ 5 x 10~ 13 along this strip. 
It would be possible to increase 6 Li/H to ~ 10~ 12 by choosing \i somewhat smaller, in which 
case the WMAP-compatible strip would be at smaller m^j. In that case, the LHC would 
enforce a stronger lower limit on mo, closer to the CMSSM focus-point strip. On the other 
hand, larger values of \i yield small values of the 6 Li abundance, and we find no increase in 
the 6 Li abundance for larger tan/3. 

The right panel of Fig. [5] displays a (/i, mo) plane in the NUHM1 for tan /3 = 20, Aq = 
and fixed m 1( / 2 = 500 GeV, at the lower end of the range allowed by the LHC and other 
data for m < 1000 GeV. In this case, there are near-vertical WMAP-compatible strips 
in transition region at ~ 300 GeV, where 6 Li/H approaches 10 -12 . There is also a 
WMAP-compatible strip near fi ~ 1000 GeV that parallels the region without a consis- 
tent electroweak vacuum (here caused by m\ < 0), where 6 Li/H is again somewhat below 
10~ 12 . Connecting these two regions is a co-annihilation segment at m ~ 100 — 200 GeV 
where the 6 Li abundance is relatively small. We have explored several other NUHM1 (/i, mo) 
planes, finding that increasing m^ decreases the attainable value of 6 Li/H. We have also ex- 
plored several other projections of the NUHM1 and NUHM2, including (m^, m^), (m^, m ), 
(/i, mA) and (mi,m2) planes, without finding values of 6 Li/H above 10 -12 . 
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m 1/2 (GeV) ,x (GeV) 

Figure 5: Left: The NUHM1 (m^, mo) plane for fi = 250 GeV, Aq = and tan /3 = 10, and 
Right: the NUHM1 (/i,m ) plane for mi/2 = 500 GeV, A = and tan (3 = 20, displaying 
contours of the 6 Li abundance including the effects of late-time XX annihilations. Contours 
of the 6 Li abundance are coloured light blue, and the WMAP- compatible strips of parameter 
space are shaded dark blue. The other shaded regions and lines have the same meanings as 
in Fig. [JJ 



4 Summary and Conclusions 

We have demonstrated in this paper that in both the CMSSM and the NUHM1 it is pos- 
sible that late neutralino LSP annihilations may enhance significantly the cosmological 6 Li 
abundance, without affecting significantly the BBN abundances of the other light element 
Deuterium, 3 He, 4 He and 7 Li [35]. This enhancement may be up to two orders of magnitude, 
yielding 6 Li/H ~ 10 -12 compared to the BBN value ~ 10" 14 . 

As we have shown, this enhancement occurs typically when the neutralino LSP is rel- 
atively light and has a large annihilation cross section, as occurs when the LSP contains 
a strong Higgsino admixture. This phenomenon appears, in particular, in the focus-point 
region of the CMSSM and in transition regions of the NUHM1. 

While interesting for the debates on the astrophysical Lithium abundances, this en- 
hancement falls short of resolving by itself the cosmological 6 Li problem. Further work 
could include a more exhaustive study of other supersymmetric models, to see whether they 
could reconcile a larger enhancement with the available theoretical, phenomenological, ex- 
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perimental and cosmological constraints. Alternatively, is it possible that the height of the 
6 Li plateau may receive contributions from other sources such as an early generation of 
stars, or might the height of the 6 Li plateau be over-estimated, perhaps because of connec- 
tive processes involving 7 Li [20]? It is clearly desirable to pin down more definitively the 
magnitude of the 6 Li problem by establishing more solidly the existence and height of the 
inferred 6 Li plateau in halo stars. However, it already seems that a substantial enhancement 
of the standard homogeneous BBN prediction for 6 Li/H might be an interesting signature 
of supersymmetric models. 

Finally, our work illustrates in detail the more general point that 6 Li production should 
play a role in-and thus probe-any WIMP dark matter scenario involving hadronic annihi- 
lation products [3T| . Specifically, we have seen that 6 Li production is essentially guaranteed 
provided there are nonthermal nucleons injected with kinetic energies > few MeV. We also 
find that the level of 6 Li abundance due to residual annihilations is model-dependent, in our 
case spanning a range from 100 times the standard yield down to an unobservable perturba- 
tion to this level. The lessons for WIMP modelers would seem to be that 6 Li observations 
already provide important constraints which one must test against, and that a confirmed 
detection of primordial 6 Li-particularly if it is above the standard level-will likely shed light 
on the details of the nonstandard physics which produced it. 
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